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We studied the feasibility of measurement of Higgs pair creation in a 77 collider. We 
found the optimum collision energy is around 270 GeV from the sensitivity study with 
Higgs boson mass of 120 GeV/c 2 . Main backgrounds are 77 — »- WW, 77 — > ZZ, and 
77 — > bbbb at the optimum collision energy. The preliminary analysis shows 77 — >• 
HH could be observed with the statistical significance ~ 5cr when we chose correct 
assignment of a track by using color singlet information. 

1 Introduction 

A Photon Linear Collider (PLC) has been 
considered as a possible option of the Inter- 
national Linear Collider (ILC). In the PLC, 
high energy photons are generated by using 
inverse Compton scattering between elec- 
trons and laser photons. An outline of PLC 
is shown in Figure 1. Details of PLC are 
found in [1]. 

In this study, we investigated the feasi- 
bility of studying the self-coupling of Higgs 
boson. The Higgs self-coupling constant A 
can be represented as A = A SM (1 + 5 k) 

which contributes via the diagram shown in Figure 2 to the Higgs pair production in photon- 
photon collision. A SM is the Higgs boson self-coupling constant which is included in the 
Standard Model, and 5k is the deviation from the Standard Model. 




Figure 1: An outline of PLC. Positron beam 
replaced to electron beam. 
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Figure 2: A diagram of Higgs pair production by 77 collision via the self-coupling. 

Theoretical studies have been performed for this process by several authors [2, 3, 4]. How- 
ever, the cross-section of Higgs self-coulping is very small, experimental feasibility against 
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huge backgrounds is yet to be studied. In this study, we investigated feasibility to observe 
Higgs pair production in high energy photon collision by Monte-Cairo simulation using a 
realistic luminosity spectrum based on a set of PLC parameters. 



2 Optimization of 77 collision energy 



In order to decide the optimum collision energy, we defined the statistical sensitivity for the 
5k as; 

\N(5k) - N SM \ L \v(t(5k) - wsm\ 



sensitivity 



^L(t]g(5k) +T1BGCTBG) 



where, N(5k) is the expected number of events as a function of 5k, and iV$M is the expected 
number of events from the Standard Model, a (5k) and <jsm are the cross-section of Higgs 
boson production as a function of 5k and for the Standard Model. L, 77, 7?bg, and (Jbg are 
the integrated luminosity, detection efficiency for signal, detection efficiency for background, 
and the cross-section of background processes, respectively. For 77 = 1 and t]bg = 0, the 
sensitivity is written as; 



sensitivity = 



a (5k) - ct S m| 
\/cr(5K) 



Figure 3 shows the sensitivity plot as a function of 
the energy of 77 interaction, v /s 77 , for the Higgs bo- 
son mass of 120 GeV/c 2 with integrated luminosity 
of 1000 fb _1 . The cross-section of signal was calcu- 
lated by using the formula given in [3] for the case of 
5k = ±1 as indicated in Figure 3. From this result, 



we set 



= 270 GeV for this study. 
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3 Beam parameters 



Figure 3: Sensitivity plot for the 
function of y/s^. Green line shows 
the optimum energy, 270 GeV. 



The parameters for the electron and laser beams are 
chosen to maximize 77 luminosity around ■ s fs 11 = 
270 GeV based on TESLA optimistic parameters [5], 
where we assumed the same parameters with the 
TESLA except for the electron beam energy (Table 

1). The luminosity distribution for the parameters was simulated by CAIN [6] and is shown 
in Figure 4. 



4 Signal and backgrounds 

Figure 5 shows the cross-section for various 77 and e + e~ processes as a function of the 
center mass energy. It indicates that the 77 — >> WW and 77 — >> ZZ processes will be the 
main backgrounds at = 270 GeV. 

Since the main decay mode of the Higgs boson of 120 GeV/c 2 is H — >• bb (branching ratio 
~ 0.68) [7], we consentrated on the case of 77 — ^ HH — ^ bbbb. Therefore, the 77 —> bbbb 
process is also considered as a background source. 
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E e [GeV] 


190 


TV/10 10 


2 


cj z [mm] 


0.35 


l^x/y [10" 6 mrad] 


2.5/0.03 


P x/y [mm] (D IP 


1.5/0.3 


a x/y [nm] 


96/4.7 


Al [nm] 


1054 


Pulse energy [J] 


10 


x — 4ouE e /ml 


3.76 




45 90 140 180 230 270 320 360 

JS^(GeV) 



Table 1: The parameters of electron and laser Figure 4: Luminosity distribution simulated 



beams. 



by CAIN. 



-|-^-h 1 ■ i ■ i i r i i f i i i ) i i ■ ] i ■ i | r i i f i i i | 

7T— *mV>>~ 




Figure 5: Cross-sections of processes as a function of center mass energy. The red line shows 
the optimum collision energy, 270 GeV. 



Then we estimated the number of events by convoluting the luminosity distribution, i.e.: 

dL 



Ch \ / S syy 



We caluculated 77 — >• HH by using the formula in [3, 4], 77 — >• WW by using HELAS [8], 
77 — >• ZZ by using gamgamZZ-code [9, 10], and 77 bbbb by using GRACE [11]. The 
numerical integration and event generation was performed by a Monte- Cairo integration and 
event generation program BASES/SPRING [12]. We expect 16 events/year for 77 — >> HH, 
1.462 x 10 T events/year for 77 — >> WW, 1.187 x 10 4 events/year for 77 —> ZZ. For 77 — >> bbbb, 
5.872 x 10 4 events/year is expected for events with bb mass of grater than 15 GeV/c 2 . 
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5 Simulation and analysis 



JSF (JLC Study Framework) [13, 14] was used as the simulation framework in this study. 
Pythia [15] was used for the parton shower and hadronization. A fast simulator [14] was used 
for the detector simulation instead of full simulation. We generated 5 x 10 4 Monte-Cairo 
events for 77 -> HH, 7.5 x 10 7 for 77 WW, 1 x 10 6 for 77 -> ZZ, and 1 x 10 6 for 
77 — >• 6666. 

First, we applied the forced 4-jet clustering to reconstruct the jets by using JADE clus- 
tering algorithm [16]. In order to select proper combination of reconstructed jets, xf s 
(i = H, W, Z, bb) were calculated as; 



mm 



(Mi -MO 2 , (M 2 -M i y 



where Mi and M2 are reconstructed particle mass. Mi (i = H, W, Z, bb) are the mass of 
Higgs boson, W boson, Z boson, and the invariant mass of bb (10 GeV). a^i = H, W, Z, bb) 
are the mass resolutions for H, W, Z , and 66, respectively. The jet combination which has 
the least x? was considered as the most probable combination. 

In order to distinguish 6 quarks, we used "nsig" 
method in this study. Figure 6 shows the outline of 
"nsig" method. For each charged track in a recon- 
struct jet, Af S ig = L/<jl was calculated, where L is 
the closest approach to the "Interaction Point" of 
the track in the plane perpendicular to the beam, 
and (Jl is its resolution. Then N Q ^ v (a), the num- 
ber of tracks which have N s \ z > a, is calculated for 
each jet. 

Before optimizing selection criteria, we applied 
the pre-selection: 




Figure 6: An outline of "nsig" method. 
6-hadron is generated at "Interaction 
Point" and decay at "Decay of 6- 
hadron". Arrows and dotted lines are 
the particle tracks and extrapolate par- 
ticle tracks. 



• A / jet(^offv(3.0) > 1) > 3, 

• A / j et(^offv(3.0) > 2) > 2, 

where A7j et (Af o ff v (3.0) > 6) is the number of jets in which Af o ff v (3.0) is greater than or equal 
to 6. 

Then we applied the Neural Network analysis to optimize the event selection criteria. 
JETNET [17] was used to train and to evaluate performance of the Neural Network. The 
Neural Network was trained for each background separately so as to maximize the statistical 
significance E defined by 

^ ^signal 

yj Signal + AbG 

where A" s i gna i and A^bg are the number of signal and background events. Table 2 shows the 
cut statistics with JADE clustering. From Table 2, the E was calculated to be 1.17<j. 

To investigate possible improvement of event selection and background suppression, we 
chose correct assignment of tracks to each jet as well as correct selection of jet pairs from 
a parent particle by using color singlet information from event generators. It was applied 
to 77 — > HH, WW, ZZ, but not to 77 —> bbbb, because color singlet combination for 66 
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HH 


WW 


zz 


bbbb 


expected events 


80 


7.31 x 10 7 


59350 


293600 


pre-selection 


47.93 


81655 


5167 


84491 


Wfilter 


12.34 


8.772 


193.4 


568.4 


bb filter 


8.238 





84.40 


13.21 


Z filter 


4.994 





7.359 


5.872 



Table 2: Cut statistics with JADE clustering. 



pair in bbbb event is not trivial. Table 3 shows the cut statistics after optimizing the Neural 
Network for perfect jet clustering. From Table 3, the E was calculated to be 4.95<r. This 
result indicates that the 77 — >> HH could be observed with ~ 5a with the integrated 
luminosity corresponds to 5-year operation of the PLC, if we could improve performance of 
the jet clustering and 6-tagging algorithm. 





HH 


WW 


ZZ 


bbbb 


expected events 


80 


7.31 x 10 7 


59350 


293600 


pre-selection 


46.64 


55836 


4172 


84491 


Whiter 


38.58 


4.873 


98.84 


2179 


bb filter 


34.50 


2.924 


27.76 


2.642 


Z filter 


33.06 


2.924 


5.935 


2.642 



Table 3: Cut statistics with perfect clustering. 



6 Summary 

We studied the feasibility of measurement of Higgs pair creation in the PLC particulary 
against large background processes such as 77 — >> WW, 77 —> ZZ, and 77 — >> bbbb. The 
optimum collsion energy was found to be 270 GeV from the sensitivity study with Higgs 
boson mass of 120 GeV/c 2 . The result is preliminary in a sense that systematic error of the 
Neural Network analysis is yet to be evaluated, but it indicates that 77 — >• HH could be 
observed with ~ 5<r with the integrated luminosity corresponds to 5-year operation of the 
PLC if we could choose correct assignment of tracks by using color singlet information. 
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